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Abstract Barchan dunes emerge as a collective phe- 
nomena involving the generation of thousands of them 
in so called barchan dune fields. By measuring the size 
and position of dunes in Moroccan barchan dune fields, 
we find that these dunes tend to distribute uniformly in 
space and follow an unique size distribution function. 
We introduce an analytical mean-field approach to show 
that this empirical size distribution emerges from the 
interplay of dune collisions and sand flux balance, the 
two simplest mechanisms for size selection. The ana- 
lytical model also predicts a scaling relation between 
the fundamental macroscopic properties characterizing 
a dune field, namely the inter-dune spacing and the first 
and second moments of the dune size distribution. 
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1 Introduction 

A first glance over an extensive desert area shows not 
only that dunes are ubiquitous and present well-selected 
shapes, but also that they typically emerge in groups 
with a very well defined characteristic dune size and 
inter-dune spacing, forming fields of up to several thou- 
sands of dunes (see Figs.QJ.-d). These observations nat- 
urally rise questions concerning the way dunes distribute 
throughout the deserts. What are the mechanisms of 
the size selection process behind such uniformity? Do 
the size distributions of such dune fields follow a simple 
unique function or do they depend on the local condi- 
tions? 

Sand dunes have been intensively studied in the 
last years. It is now well-understood what are the fun- 
damental laws underlying the emergence of one single 
barchan dune and what mechanisms maintain its shape 
while moving 1,2 . For instance, barchans occur in ar- 
eas with unidirectional wind and low sand availability. 
The influence of the geographical constrains and the 
external physical conditions [2|[3] , of the dune-dune in- 
teractions [HEIE] and even of the emergence of vege- 
tation covers [7] in the dynamics and morphology of 
single dunes were quite well-established with the help 
of dune models [TJ1H1E]. There are also a few studies of 
entire dune fields [ MTUlfTTlll^] , but a simple theoretical 
understanding of the size selection process within dune 
fields has still not been achieved. 

In this work, we present a first answer to this prob- 
lem. First, we show that, while a single dune is suitably 
characterized by its width w [2lll2j , an entire dune field 
contains dunes with different sizes following a unique 
distribution (Fig. [Tf-g) . Consequently, the correspond- 
ing average (w) and standard deviation S = y (w 2 ) — (w) 7 
are suitable properties to characterize the field. Ad- 
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Fig. 1 (Color online) (a-d) Details of barchan dune fields in Morocco, Western Sahara. The number of measured dunes is 1295 (a), 
1113 (b), 1947 (c) and 1630 (d), covering areas of ~ 3,7,12 and 60 km 2 and with average dune sizes of 17m, 27m, 42m and 86m 
respectively. The size of a barchan dune is characterized by its width w (e). In all pictures, the North points up. Images provided by 
GoogleEarth. (f ) Probability density function (PDF) of the dune size for the measured Moroccan dune fields (symbols) and the best 
fit using the analytical solution (solid lines) given by Eq. {T}. (g) Cumulative distribution function (CDF) with the analytical solution 
(solid lines) and the log- normal straight-line for reference (dashed lines). The relative broadness S/(w) is given by the inverse of the 
slope of the CDF. After rescaling the dune sizes as (w/{w)Y with t = 2.9%/{w)/(L), all PDFs and CDFs in (h) and (i) respectively 
collapse, uncovering a scale invariance between the size distributions of different barchan dune fields. In (h) the analytical distribution 
(see text) is also shown (dashed line), as an eye-guide. 



ditionally to these two quantities, we show that the 
inter-dune spacing L is also a property with character- 
istic values (L) and therefore also able to characterize 
the field. Second, using numerical simulations, we show 
that collisions between dunes play a crucial role in the 
selection of a characteristic dune size. Finally, from a 
mean-field approach that couples the effect of dune col- 
lision with that of sand flux balance, we derive the size 
distribution function (shown in Fig. Hf-g) and a scaling 
relation between the three properties of the field, (L), 
(w) and S. 

2 Measurements of the dune size and 
inter-dune spacing distribution 

We start by measuring the width (Fig[T^) of more than 
5000 dunes composing four dune fields located in the 
Western Sahara (Figs. QjL-d). For all fields, the dune 
width distribution exhibits a unique function, apart 
small deviations at the extremes (w < fOm and w > 
0.5km[T2]), as shown in FigQJ-g. This function will be 
derived later. The relative broadness S/(w) scales with 



the relative inter-dune spacing as yj (L) / (w), as shown 
in Figs. [TJi-i. This scaling law relates the spatial dis- 
tribution of dunes and their size distribution and will 
also be deduced later as a result of the size selection 
model we propose. Therefore, the mechanisms leading 
to such distributions should not depend on the absolute 
size and inter-dune spacing of the dunes involved. In- 
stead, they should depend on the relative dune size and 
spacing. In other words, they should be scale invariant. 

In barchan dune fields, the sand flux balance on sin- 
gle dunes leads to an instability in the dune size [5J. 
Dunes nucleate as a consequence of the sand accumula- 
tion along the field and no characteristic size emerges [6j 
10 . However, the peaked size distributions in Fig. QJ are 
found in dune fields where sand flux balance is not the 
only process mediating the size of dunes. Since barchan 
dunes move over the field with velocities that strongly 
depend on the size (v oc 1/w) [1] , collisions are ubiqui- 
tous in such fields, turning out to be another relevant 
process for dune size alterations [Blll2j. Therefore, the 
dune size distribution should be determined by the com- 
petition between the balance of sand flux on a single 
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Fig. 2 (a) The PDF of the inter-dune spacing clearly shows a characteristic value (L) for each dune field depicted in FigUfa-d). The 
inter-dune spacing L around a given dune is defined as the square root of the empty area of a polygon formed by the centers of the 
nearest dunes, one in each quadrant of the Cartesian frame reference centered at the dune, (b) PDF of L as a function of the dune 
width w for the first dune field a. From the contour plot, one defines the characteristic inter-dune spacing (L) (solid line), taken as the 
average over the highest frequency region, which is independent of w (see text), (c) After rescaling L by the mean dune size (w), not 
all PDFs peak in the same relative inter-dune spacing (L)/(w). (d) However, the curves collapse after rescaling L by the expression 
S 3 /(w) 2 , where S represents the standard deviation of the size distribution (see text). 



dune and the collisions between neighboring dunes [5J 
151110]. Recently a third size selection mechanism was 
discovered, that involves the calving of large dunes due 
to wind fluctuations [12,13J. This is a complex scale 
dependent process, relevant for fields with large dunes 
and should not be responsible for the distributions here 
addressed. Therefore, we will not consider it. 

Two important aspects must now be addressed to 
proceed further. First, we notice that in the absence of 
motion, i.e. in static fields where no collisions can occur, 
one dune grows only if its neighboring dunes shrink [14 , 
due to sand flux balance and mass conservation. Con- 
sequently, if collisions do not occur the inter-dune spac- 
ing L between neighboring dunes would scale with the 
dune size. Our empirical data however, shows a rather 
different behavior. For each dune field in Figs.QJi-d, the 
inter-dune spacing L between each dune and its neigh- 
bors as a function of the width w, distributes parallel 
to the w-axis (Fig.[2b>). Thus, contrary to the situation 
without collisions, here the inter-dune spacing can be 
taken as its characteristic value, say (L), as shown in 
Fig. [2^ and [2Jd. This empirical result is a clear sign of 



a richer internal dynamics in dune fields where colli- 
sions play an important role. Indeed, due to collisions, 
small dunes are continuously emerging from larger ones 
[121113] , destroying any simple correlation between dune 
size and inter-dune spacing, and therefore leading to the 
observed spatial uniformity. 

The second aspect is that (L) does not follow a sim- 
ple scaling with the average dune size (w) of the field 
(Fig. [2b) but a more complex one (as S 3 /(w} 2 , Fig.[2U) 
which also involves the standard deviation S of the size 
distribution. As will be shown later, this scaling nat- 
urally emerges from the coupling of binary dune colli- 
sions and sand flux balance in the size selection process. 

3 Binary collision dynamics 

With these two empirical findings we proceed by study- 
ing the effect of collisions alone in the size distribution 
of a dune field. Using an established dune model [IBS], 
we simulate ideal binary collisions, under open bound- 
ary conditions and constant wind (Fig[5feo), extracting 
a simple collision rule, i.e. a phenomenological function 
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that relates the relative dune size rt after the colli- 
sion between two dunes, and the corresponding initial 
relative size and offset 9i [15] , In accordance to re- 
cent underwater experiments [6] and our simulations, 
for most initial conditions the number of dunes is con- 
served and the total sand volume change is negligible. 
Furthermore, as seen in Fig. [3Jd, in most of the cases 
the collision increases the relative dune size (ry > r^), 
redistributing sand from large to small dunes, in sharp 
contrast to the flux balance, which accumulates sand 
on large dunes [15]. 
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Fig. 3 (a) Different outcomes of simulated binary collisions, 
from coalescence (when both dunes merge) on top, to a situation 
where the volume of the smaller dune increases (decreases) after 
the collision (middle and bottom sections, respectively), (b) Col- 
lision rule for binary collisions that conserve the number of dunes. 
Dots represent numerical simulations and dashed lines the cor- 
responding surface fit rj(ri,9i). The curve rj(r{,9i) = (solid 
line) separates two regimes: one with rj > rj, where collisions 
redistribute sand and another with rt < ri due to accumulation 
of sand. 



With this collision rule we study the influence of 
collisions on the dune size distribution neglecting sand 
exchange through sand flux between them. A large sam- 
pling statistics of about 10 000 dunes is considered, 
where we assume that each pair of dunes collides with 
the same fixed probability following the collision rule 



introduced above (Fig. [3)3 ) ■ This assumption does not 
consider the spatial distribution of dunes and thus we 
neglect the spatial correlations between dune sizes and 
positions. As a result of binary collisions, arbitrary ini- 
tial size distributions converge to a stationary Gaussian 
distribution P co i(w), illustrated in Fig.0Ji. We find that 
the relative standard deviation <j co ij (w) co i of such dis- 
tribution, where {w) co i is the mean size, is constant for 
different initial conditions (Fig. [4})) and its value is fixed 
by the parameters of the dune model, which were chosen 
to reproduce the morphology of Moroccan dunes 0[T6j , 
which in turn determine the dune morphology and thus 
the phenomenological collision rule [15] . The constant 
value of the relative standard deviation is consequence 
of the scale invariance of our collision dynamics that de- 
pends on the relative dune size instead of the absolute 
size. On the contrary, with changing wind conditions, 
the dune scale may influence the collision outcome and 
also several dunes may emerge from a single collision 

Until now, two preliminary conclusions can be stated. 
First, the measured barchan dunes are approximately 
uniformly placed over the deserts, with characteristic 
values of L, and the different dune fields follow a com- 
mon size distribution with a simple scaling with the 
mean dune size. Second, neither sand flux alone does 
contribute to the size selection [6lfT0] nor collisions alone 
can be responsible for the skewed measured non-Gaussian 
distributions (FigQJ. What results from the interplay 
between both mechanisms will be now derived from a 
mean-field approach, leading to a master equation for 
the size distribution P(w, t) and to a simple scaling be- 
tween all the three properties, (w), S and (L). 

4 Master equation 

When only ideal collisions occur, Fig. [4J: shows that 
both, the mean (w) and the standard deviation a of 
the size distribution exponentially relax toward their 
respective equilibrium values, (w) co i and a co i. By using 
the definition (w)(t) = J dw w P{w 1 i) and the expo- 
nential relaxation d(w}/dt — ((w) co i — (w))/t c , the size 
distribution P(w, t) obeys in first approximation the 
dynamical equation dP/dt = {P co i — P)/t C i where t c is 
the characteristic relaxation time (in units of number 
of collisions per dune) of P(w) towards the equilibrium 
Gaussian distribution P co i when only collisions are con- 
sidered. 

When both collisions and sand flux balance are con- 
sidered, the total temporal derivative of P(w, t) has now 
two separated terms, the partial temporal derivative 
dP/dt and the term arising from the volume change 
rate V = wdV/dw, due to flux balance: VdP/dV = 
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Fig. 4 (a) Snap-shots of the evolution due to binary collisions of a uniform dune size distribution (dashed-lines) towards a Gaussian 
(solid line), (b) Linear relation between the two first moments of the observed Gaussian, for different initial volumes of sand available, 
yielding a constant relative standard deviation cr co ;/{ui} ( . oi . (c) The first two moments (w) (squares) and cr (circles), exponentially 
relax toward their equilibrium values (w) ca i and u co i (see text). Time units are number of collisions per dune. 



wdP/dw. Assuming the existence of a steady state and 
using the empirical fact that V oc w 3 and V oc Qw, with 
Q denoting the saturated sand flux over a flat bed [12] , 
the master-equation yields 



sonable approximated as 



t c dP(w) 
t s dw 



(w) 



(P col (w) - P(w)) , 



(1) 



where there are three parameters, namely the relative 
deviation o C oi/{w)coi of P co h determined by the colli- 
sion model, the characteristic size {w) co i and the ratio 
t c /t s . Time t s is the characteristic time associated to 
the change rate of the dune size due to the sand flux 
balance, defined as t s — a(w)1 ol /Q with a as a con- 
stant. From Eq. (fj]) one concludes that when collisions 
dominate in the selection of dune sizes, t c <C t s and 
consequently the distribution converges to the Gaus- 
sian P co i. When the opposite occurs, with the sand flux 
balance being the relevant process, P(w) deviates from 
P co i(w) the more the larger t c /t s is. 

The solution of Eq. (fT]) is plotted in Fig. QJ and [Tg 
(solid lines) for each dune field, with (u>) co ; and t c /t s as 
fit parameters. The value of a co i is taken from Fig. 2jo. 
Apart extreme points, the solution fits well the empir- 
ical distributions, with first and second moments rea- 



(w) 
S/(w) 



0.8(w)coi(VtJts + l) 

Qm^/tjr s /(^/i~jr s + i) 



(2a) 
(2b) 



in the range t s < 5t c . From Eqs. one sees that 
the characteristic size (w) co i determined by the colli- 
sions dynamics is in fact the only characteristic size in 
the system. Moreover, Eq. (|2b[) shows that the relative 
standard deviation S/ {w) is given by the ratio t c /t s 
and thus describes a measure for the competition be- 
tween sand flux balance (t s ) and collision (t c ) processes 
for the dune size selection. For instance, dune field in 
Fig- Ufa has a large value t c /t s = 4.3, indicating that the 
sand flux balance is the most important size selection 
process, while the dune field in Fig. QJ: has t c /t s = 1.7 
indicating more relevance from collision processes. 

Furthermore, the ratio t c /t s is not an independent 
parameter since t c must be proportional to the col- 
lision time t co i, defined as the average time for two 
dunes to collide. This collision time is determined as 
the quotient between the inter-dune spacing (L) and 
the average relative velocity between two dunes, (v) r oc 
J °° dw\ dw2P{wi)P(w2)(v(wi) —v(w2))- Since the 
dune velocity follows v oc Q/w and, within some wide 
range of sizes, P(w) can be well approximated by a log- 
normal distribution (see Figs. QJ and[Tg), one obtains 
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(v) r cx (v}(T£ with (v) cx Q/(w) the dune average ve- 
locity and (T£ the standard deviation of the log-normal 
distribution (adimensional), yielding t c cx ■ Sub- 

stituting t c and Eq. (|2a|) into Eq. (|2b|) and taking the 
first-order approximation cr^ ~ S/(w), we arrive at 

(S/M) 3 ^ a <L}/<». (3) 

As shown in Figs. [T]i-i and[2tlj where one empirically 
obtains e = (^/(w}) -1 = 2.9{/{w)/{L), the scaling in 
Eq. ([3]) indeed describes the measurements, with the 
constant a = 2.9 -3 = 0.041, independent of the model 
parameters. 

5 Conclusions 

In summary, due to the dynamical nature of barchan 
dunes, the size distribution is intrinsically linked to the 
spatial distribution in such a way that sparse fields 
have a broader size distribution, while dense ones have 
narrow distributions. We have shown that the relative 
dune size distributions of Moroccan dune fields col- 
lapse into an unique distribution function, and that 
dunes are uniformly distributed with a characteristic 
inter-dune spacing that obeys a simple scaling law. By 
using a master-equation approach with a simple colli- 
sion rule, we showed that the simplest processes behind 
the change of the dune size occurring in dune fields, 
namely ideal binary collisions and flux balance, are able 
to properly determine the size distribution function. 
Which mechanisms are behind the local selection of 
the specific size scale of a dune field remains an open 
question, since it involves not only binary collisions and 
flux balance under a stationary wind, but detailed pro- 
cesses in real changing wind conditions, i.e. calving [12[ 
fl"3"] . that can locally change the dune size and are not 
included in the analysis we have presented. 
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